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QUANTITATIVE  INTERPRETATION  OF  DISLOCATION 
CONTRAST  IN  ELECTRON  MICROSCOPY 


By 


Hans-Jurgen  Kestenbach 


Major  Department:  Metallurgical  and  Materials  Engineering 

Electron  microscope  images  of  dislocations  in  aluminum  have 
been  studied  by  quantitative  comparison  with  calculated  intensity 
profiles  predicted  by  dynamical  electron  diffraction  theory. 


could  be  derived  a. 


true  dislocation  line  position  agreed  with  calculations.  Rather  strong 

both  systematic  and  weak  nonsystematic  reflections.  Image  sensitivity 
to  several  contrast  parameters  has  been  investigated,  and  it  was 

contrast  is  a function  of  foil  thickness,  in  agreement  with  measurements 
from  the  Intensity  distribution  aeross  thickness  extinction  fringes. 


For  Burgers  vector  determination,  a g 
which  provided  immediate  information 
during  observation,  and  which  can  als 


t the  character  of  a dislocation 


When  a crystalline  specimen  is  imaged  is 
. electron  microscope,  contrast  is  produced  mainly  by  electron  diffrac- 
tion. thus,  local  regions  of  strain  can  be  distinguished  from  the 
of  perfect  crystal,  since  they  will  either  enhance 
:he  diffracted  intensity.  It  is  said  th8t  such  local  dlsplace- 

is  realized  when  the  first  electron  microscope  Images  of 

obtained  in  1956,  Independently  by  Hirsch,  Horne  and 
Whelan  ^ and  Bollmann. ^ Since  then,  strain  contrast  has  played  the 
major  role  in  electron  microscopy  as  a tool  for  metallurgical  research. 

Within  a few  years,  the  contrast  mechanism  had  been 
theoretically  established,  by  Whelan  and  Hirsch^  and  Howie  and 

the  more  sophisticated  dynamical  theory  of  elec 


identify  Burgers  vectors  of  both  t 


nt  precipitates.^  For  an 
le  reader  is  referred  to  the  book  by 

it  led  to  the  possibility  of 

loops  of  vacancy 
considerable 


tn  expended  ii 


for  the  study  of  radiation  damage. 

Despite  the  considerable  effort  which  hi 

sophisticated  and  quantitative  contrast  theory  are  use< 
interpret  experimental  data  only  in  a qualitative  team 
Burgers  vectors  are  identified  when  dislocation  images 
under  particular  imaging  conditions,  and  strain  fields  of  small 
defects  are  analyzed  by  determining  the  crystallographic  directions  of 
lines  of  symmetry  in  their  images.  Yet,  the  dynamical  contrast  theory, 
at  least  in  principle,  is  capable  of  predicting  the  intensity  level  at 

has  been  chosen.  It  should  then  be  possible,  in  suitable  applications, 

which  allows  accurate  and  quantitative  matching  between  experimental  am 
theoretical  information.  This  constitutes  the  main  objective  of  the 


The  study  of  contrast  from  inclusions  by  Ashby  and  Brown'  ' ‘ 1 
first  attempt  toward  a more  quantitative  evaluation  of 
strain  field  Images,  with  the  aim  of  measuring  the  magnitude  of 
coherency  strains  for  precipitates  of  various  shapes.  The  authors 
calculated  image  width  at  a given  intensity  level  (with  respect  to 
background)  as  a function  of  strain,  and  compared  experimental  Images 
to  intensity  scales  which  were  prepared  on  photographic  places  in 
order  to  measure  the  width  at  a particular  intensity  level.  A more 


lattice  defects  whicl 


een  the  image  simulation  technique,  originally 
n directly  be  compared  with  experimental 


available,  this  aetl 
Specifically,  it  dot 
of  position  which  wc 


Alchough  utilizing  the  t 

i on  qualitative  interpretation. 

.f  quantitative  predictions  about 
the  strain  field  are  to  be  verified.  The  research  work  to  be  presented 

calculating  one-dimensional  intensity  profiles,  and  has  been  directed, 
whenever  possible,  toward  the  substitution  of  experimental  intensity 


Intensity  profiles  ai 


photographs  by 


to  be  more  reliable,  since  the  variable 
processing  are  no  longer  involved.  In 
successful,  can  be  extended  to  Include 
Investigate  the  effects  of  inelastic  sc 

s quite  plausible  and  also 


s expected 


addition,  this  second  method,  if 
ilectron  energy  analysis  to 
:atterlng.  An  intensity  scanning 


is  only  this  dependence  which  stimulates  metallurgical  incerest  In 
contrast  theory.  However,  calculations  have  shown  chat  strain  field 
images  are  also  quite  sensitive  to  several  ocher  variables  which  are 

position  of  the  defect  in  the  foil,  the  operating  reflection  (or 
reflections  in  many-beam  contrast) , the  precise  foil  orientation  with 
respect  to  the  electron  beam,  and  finally  two  parameters  which  describe 

contrasc  parameters  are  not  precisely  known  or  difficult  to  measure 
in  a practical  situation  (with  the  exception  of  Che  operating  reflec- 
. tion),  but  their  values  must  be  assessed,  at  least  to  some  extent, 

The  final  aim  was  of  course  to  find  imaging  conditions  under  which  all 
parameters  are  sufficiently  well  known  to  an  extent  which  would  make 
the  image  representative  of  the  displacement  field  in  Che  foil. 


la  often  difficult  in  practice  Co  achieve  optimal  cwo- 
It  is  therefore  desirable  to  know  under  what  experimental  conditions 

beam  contrast  has  been  Included  in  the  present  Investigation  to  help 
answer  this  question. 

Since  the  research  objective  as  outlined  above  consisted  mainly 
of  exploring  the  feasibility  and  accuracy  of  quantitative  image 

a known  displacement  field  could  accurately  be  recognized  by  its 


for  such  a task  because  the  lattice  distortion  due  Co  dislocations  can 
be  calculated  with  confidence  from  elasticity  theory.  The  use  of 
isotropic  elastic  theory  to  calculate  displacement  fields  led  to  the 

are,  however,  other  advantages  which  aluminum  offers  in  a study  of 
diffraction  contrast.  First,  its  low  scattering  power  for  electrons 
allows  the  observation  of  dislocations  in  rather  chick  foils  in  which 
surface  effects  on  the  displacement  field  will  play  a lesser  role. 

more  realistic  in  crystals  containing  light  atoms. 

complicated  contrast  mechanism  can  yield  quantitative  information  about 


profiles  of  di; 


e accuracy  of  quantitative  image  Interpretation, 
determining  Che  exact  shape  of  dislocation  profiles,  including  many- 
The  following  chapter  will  present  an  outline  of  Che  dynamical 
which  dislocation  contrast  has  been  calculated.  Details  relating  to 
describes  the  intensity  scanning  system,  specimen  preparation,  and  the 
parameters.  Chapters  1 and  5 finally  are  reserved  for  presentation 
been  placed  on  the  formulation  of  final  conclusions.  The  appendix 


CONTRAST 


2.1  Introduction 

scattering  events,  between  high-energy  electrons  end  matter,  any- 
one o£  which  can  he  a potential  source  for  image  contrast  in  the 

electron  retains  its  original  energy,  or  is  celled  inelastic  if  the 

crystalline  specimens,  the  elastic  scattering  is  governed  by  Bragg's 

plays  the  dominant  role  in  the  formation  of  image  contrast.  There- 

the  past  are  based  on  the  theory  of  electron  diffraction. 

diffraction,  the  situation  in  the  electron  case  is  quite  different 
from  a practical  point  of  view,  the  reason  being  the  much  higher  rate 


Bragg  orientation,  * 


diffracted  beams  Mill  be  present. 

.s  traversing  t) 

The  Intensity  carried  in  any  particular  beam  direction  will  be  • 

the  problem  is  to  calculate  che  intensity  of  the  did 
of  position  in  che  foil. 

the  historical  development  of  which  has  been  briefly  outlined  in 

There  are  two  basically  different  approaches  for  deriving 

The  first  approach  has  been  applied  by  Bowie  and  Whelan^  co  the  two- 
beam  theory,  and  by  Howie  and  Basinski,  based  on  a method  devised 
by  Takagi,^*18^  to  the  many-beam  theory.  It  has  been  called  the 


am  (bright 


c formation, 
if  diffraction, 


INCIDENT 

ELECTRONS 


Schematic  representation  of  diffraction  contrast 
produced  by  a dislocation.  A dark  field  image 

Intensity  is  indicated  by  arrow  length. 


Che  formal  equivalence  of  Che  two  approaches  in  Che  cwo-beam  cheory. 

In  che  following  section,  Che  oucline  of  the  Bloch  wove  treat- 
ment is  briefly  discussed,  which  is  more  involved  mathematically , but 

addition,  ic  is  shown  how  ics  equations,  in  the  general  many-beam  case, 

reasons.  There  is  on  important  difference  becween  these  two 
categories  which  should  be  emphasized  at  this  point.  Mathematical 

a physical  process,  and  their  validity  can  be  assessed  on  a theoretical 
basis  without  recourse  to  an  experiment.  On  the  other  hand,  physical 
assumptions  necessarily  change  any  theory  from  a rigorous  treatment 


The  most  satisfactory  and  rigorous  treatment  of  che  dynamical 


72«?)+S5JS[e  + v<J>]*<?>-0 


j allowed 


which  have  wave  vectors  St,  corresponding  to  the  energy  levels 


The  crystal  potential  V(r)  is  usually  written  as  a three- 


where  the  U and  U'  are  the  Fourier  coefficients  related  to  the  crystal 
potential  for  elastic  and  inelastic  scattering,  respectively,  and  the 

cients  are  related  to  atonic  scattering  amplitudes  via  the  well-known 
structure  factors,  F : 


with  V representing  the  volume  of  a unit  cell  in  the  crystal.  T1 

is  more  complex  since  the  phenomenon  of  inelastic  scattering  is  ni 
yet  fully  understood  on  a theoretical  has 
in  several  materials  have  recently  been  calculated,' 
accuracy  has  not  yet  been  established.  This  point  w 

scattered  electrons  from  the  beam.  This  is  equivalent  to  the  absorption 
of  electrons,  and  U'  is  therefore  called  the  absorption  coefficient. 


orward  direction,  and  many  i 


le  inelastic  scattering  proceeds  it 


ten  incorporated  into  the  dynamical  theory. 
Substitution  of  (2.3)  and  (2.4)  into  (2.1] 


from  which  K can  be  identified  as  the  wave  vector  of  the  electrons  in 
the  crystal  after  refraction  by  Che  mean  lattice  potential,  U . 
Equation  (2.6)  represents  a set  of  n equations,  each  for  a given 
reciprocal  lattice  vector,  g,  and  containing  a sum  taken  over  all 
reciprocal  lattice  vectors,  h.  The  system  is  linear  in  its  unknown 
variables,  C.  , which  depend  upon  the  eigenvalues,  representing  those 
values  of  k for  which  the  system  has  a solution.  Each  eigenvalue,  k, , 


detail  in  Appendix  A: 


e general  solution  t 
p is  not  trivial,  it 


#<*>  * I »,  I C (k  )exp[2rl(k 
} i S S J 1 


Although  kj  and  C <k.)  are  determined  b; 
must  be  determined  from  boundary  conditions.  This 


6),  equation  (2.3)  now 
.1  depend  on  the  particular 


The  boundary  conditions  are  illustreted  in  Figure  1 where  the 
physical  situation  is  described  by  the  electron  distribution  before 

the  intensity  is  carried  in  the  incident  beam.  Equation  (2.8)  is 
then  expected  to  explain  the  situation  after  the  electrons  have  left 


the  foil.  In  both  cases,  the  electron  wave  function  of  (2.8)  must  be 
expressed  in  terms  of  plane  waves  travelling  in  certain  directions 
(K  + g)  which  represent  the  direct  beam  entering  the  foil  and  the  direct 
beam  and  several  diffracted  beams  leaving  it* 


*(r)  - l dgexp[2vi(K 


g)  *r) 


with  amplitude  6 obey  Schrodinger's  equation  in  vacuum.  The  only 
requirement  is  that  equation  (2.6)  Inside  the  foil  join  smoothly  wltl 
the  wave  function  in  vacuum  on  both  foil  surfaces.  There  is  thus  no 
assumption  required  for  the  electron  distribution  in  the  crystal,  ant 
it  is  at  this  point  that  the  Bloch  wave  treatment  is  more  rigorous 


Equating  (2.8)  and  (2.9),  the  bean  amplitudes  can  be  derived 
as  a function  of  position,  r,  which  finally  determine  the  intensity 

discussed.  The  mathematical  formulation  is  based  on  an  analytical 


Che  solution  to  equation  (2.6)  in  the  so-c. 
representation  which  derives  its  name  from 
of  Che  electron  beam  impinging  on  the  foil 

dispersion  is  governed  by  equation  (2.6),  . 


using  matrix  notation, 
led  dispersion  surface 

energy  remaining  always 
b solution  is  referred  to 


is  illustrated  in  Figure  2,  which  has  been  drawn  t 
case  of  five  systematic  reflections  being  excited 

referred  to  the  literature. 

Figure  2 shows  the  five  pertinent  branches  of  the  dispersion 
surfaces  together  with  the  five  wave  points  corresponding  to  the  fi\ 


a row  of  reciprocal  lattice  points  c 


Schematic  five-beam  case  ii 
representation.  All  circlt 

kinetic  energy  associated  ; 


dispersion  surface 


line  which  Is  perpendicular  Co  Che  foil  surface,  this  Insures  Che 
cangenclal  continuity  of  che  eleccron  wave  vector  afcer  entering  che 
cryscal.  Thus,  boundary  conditions  have  already  partially  been 
fulfilled  in  drawing  Figure  2,  K,  as  defined  by  (2.7) , is  Che  wave 

che  foil  surface),  and  represencs  che  radius  of  che  Ewald  sphere  in 

connecting  che  Bloch  wove  poinCs,  and  in  che  figure  is  almost  coinciding 
wich  che  wave  point  of  Bloch  wave  1.  y gives  the  distance  of  che  wave 

Figure  2,  a convenlenc  relationship  between  che  wave  veccors  SC  and  it. 


K2  - I®  + Jl2  % K2  - [K  + Or  - sg)cos  6g]2  % 

% 2K(sg  - y)cos  8g 


e approximate  idendcy 


practice  all  Bragg  angles  in 


be  close  Co  unity  and  can  be  neglected. 

We  are  now  in  a good  position  to  start  che  matrix  formulation 


M©-»© 


■’srsfc; 


" Vf>  Vf  but  Ia  ei 


of  (2.9)  to  equation  (2.8)  which  represents  the  total  wave  function  of 
the  electrons  In  the  foil.  If  (2.8)  is  regrouped  into  terms  correspond- 

are  found  to  be,  using  (2.10): 


■ l *,C  (y,)cxp(2riy  y) 


em  of  figure  1.  In  deriving  (2.14), 
from  Figure  2 that  K and  kj  only  differ  in  their 
difference  given  by  y, . Equation  (2.14)  can  be 


(*>-  [£]  { «P«viYjy)}  <♦) 


i represents  a 


is  a column  vi 
denote  a diagonal  m: 

Equation  (2. 


■ 


th  elements  5 ,4  ,...,«  , 

) arranged  in  the  order  j ■ 0,1,2 and 

*.  and  exp(2riy^y)  will  be  the  Jth  diagonal 
can  now  be‘  used  to  relate  beam  amplitudes  at 
parameters  ^ of  the  Bloch  waves 


<*>  • 


be  substituted 


foil  surface,  only  Che  fir3E  element  $0  in 
ic  follows  Chat  ijij  - C0(Yj).  Equation  (2.16)  c 
into  (2.15)  Co  obtain  a relationship  between  th 


crystal  slab  of  thickness 


(•>„•[-]  (•>.. 

[']  ■ H {"•«'»>*>} [«]- 


. if  a perfect  crystal  is  imaged 
directly  gives  the  amplitude  ar 
bottom  surface  in  any  desired  l: 
No  physical  assumption 
of  perfect  crystal 


ist.  However,  two  assumptions  must  be 
:y  is  presently  extended  to  apply  to  Imperfect 

crystal  potential  to  be  fixed  in  space  with 
: Che  atom  is  displaced  by  a small  distance 
so  is  the  particular  value  of  the  potential  originally  found  at 
'he  expression  for  the  crystal  potential  which  mus 
Schrodinger's  equation  (2.1)  is  Chen  given  by 


Uh,  so  is . 


M') exp [2*1  g'(r  - R<r)>] 


Substituting  (2.3)  and  (2.19)  into  (2.1),  a new  eigenvalue 
equation  follows  which  is  identical  in  fore  to  (2.11) 


M <=•> 


, - <,) 


defined  as  before,  the  prltte  indicating  that  their 
values  or  elements  will  in  general  be  different  from  those  of  the 
perfect  crystal.  [a'J  has  diagonal  elements  which  are  identical  to 


exp [-2*1 (g  - K) •$(?)) 


This  suggests  that  the  imperfect  crystal  should  be  divided  first  into 
narrow  vertical  columns  (the  column  approximation  used  by  Howie  and 

further  divide  any  column  into  small  slabs  of  thickness  oy  which  can  be 
assumed  to  represent  slabs  of  perfect  crystal,  and  which  therefore 
should  obey  equation  (2.17).  This  constitutes  the  second  assumption 


elements  of  the  diagonal  m. 


Taylor  expansion 


exp<2siyjSy)  - 1 + 2rlY^6y 


With  [c]  jyj}  [c]"1  - [ a]  from  equation  (2.11),  and  in  Che  limit 
of  infinitely  small  dy,  (2.23)  can  be  written  as  a differencial 


: (♦(?>)  - 2ri  [a(?)]  <«(y)> 
notl“  th“  57  (♦*)•* 


iow  expressed  ii 


13)  a: 


h represents 

scattering  over  a distance  dy  in  the  f 
of  quantities  which  no  longer  include 

done  numerically),  the  solutions  to  the  dispersion  equation 
to  the  general  eigenvalue  equation  (2.11)  or  (2.20)  are  no  longer 

overlooked  due  to  Che  simplification  of  the  matrix  notation  becaui 
appears  to  be  the  most  ingenious  and  puzzling  feature  of  the  Bloc 
wave  treatment  developed  by  Howie  and  Whelan.  It  significantly 
reduces  the  amount  of  computations  necessary  for  Intensity 


“1 1 Is®-******®  * w 


In  order  Co  solve  a system  of  linear  differencial  equations, 
ac  lease  one  independenc  boundary  condidon  muse  be  specified  for 
every  equacion,  which  In  chis  case  is  found  from  Che  elecCron  dlscrl- 
bucion  on  Che  encrance  surface  of  Che  foil,  and  which  is  equivalenc 
Co  Che  boundary  condidon  discussed  in  Che  derlvacion  for  the  perfect 
crystal.  At  the  foil  entrance,  all  Che  intensity  is  carried  in  Che 
direct  incident  beam  and  is  chosen  to  be  unity  for  convenience.  The 
sec  of  boundary  conditions,  equivalent  to  equation  (2.9),  thus  becomes, 
in  the  coordinate  system  of  Figure  1, 

(2.30d) 


Since  equations  (2.30a)  and  (2.30b)  in  practice  can  neither 
be  summed  over  an  infinite  number  of  reciprocal  lattice  vectors  h, 
nor  contain  an  infinite  number  of  equations  (one  for  each  reciprocal 
lattice  vector  g),  a decision  must  be  made  as  to  which  reflections 
can  be  neglected.  In  a contrast  experiment,  this  can  be  judged  from 


relative  intensity  of  diffraction 


liti: 


ia 


10X.<9)  Recently,  Metherell  and  Fisher^^  investigated  the  theoretical 
sensitivity  of  many-beem  thickness  and  stacking-fault  fringes  to 

Itao-beam  calculations  are  always  carried  out  dimensionless 

across  Che  dislocation  will  be  measured  in  angstroms,  which  means  that 
also  for  two-beam  contrast  numerical  values  become  important.  For  best 

have  therefore  been  used  in  all  two-beam  calculations.  In  many  practical 
situations,  however,  the  error  involved  in  the  measurement  of  geometrical 
be  larger  than  the  uncertainty  with  respect  to  the 

It  is  only  for  reasons  of  notation  that  the  term  € (correspond- 


le  dynamical  theory  distinguishes  bi 


phenomenon  exists  in  electron  diffraction  or  not  is  quite  immaterial 
for  the  shape  of  an  intensity  profile  which  merely  depends  on  a 
parameter  s'  representing  the  effect  of  inelastic  scattering.  The 
absorption  parameters  are  usually  expressed  as  ratios  £ /t'.  and  the 
following  definition  will  be  used  from  here  on  for  convenience: 


The  mean  absorption 


(2.27)  by  a substitution  similar  to  those  of  (2.26)  oi 
therefore  only  changes  the  overall  Intensity  level  of 


strain  field.  The  value  of  0.10  has  thus  been  assumed  quite  arbitrarily 
for  A,  expect  in  some  instances  where  profiles  have  been  calculated 
for  B > 0.10  since  the  case  A < B leads  to  the  physically  Impossible 


of  intensity  fringes  in  thick  foils.  Recently,  a 


e reduced  visibility 
d in  particular  with 


ay-beam  calculations  to  prepare  for  high-voltage 


calculated  by  Humphreys 
scattering  into  account,  whereas  t 
electron  excitation  ai 


Anomalous  absorption 

Hirsch^2^  and  by  Rad*,  2t*^  derived  fi 
electron  scattering  theory.  The  firs! 


, considering  all  three  principal 


above.  Experimental  data  from  thickness  fringe  mi 
added  as  far  as  they  were  available.  A very  careful  investigation  of 
experimental  measurements  of  absorption  parameters  has  been  undertaken 
by  Bell  ^2^  whose  data  may  be  most  trustworthy  in  Table  1. 

has  been  assumed  in  matching  calculations  for  all  reflections  encounterei 
This  was  in  part  justified  by  preliminary  computations  which  showed 
some  dislocation  profiles  being  little  affected  with  respect  to 
their  shape  by  changing  B over  a wide  range  of  values.  In  addition, 
fisher^22^  reported  that  the  function  of  B with  respect 

fringe  profiles.  However,  further  calculations  during  the  course 


ABSORPTION  PARAMETERS,  B,  FOR  ALUMINUM 


quite  influential  in  particular  contrast  situations,  sene  of  which  will 
be  discussed  in  decall  in  Section  4.2. 


in  y-direction.  The  Ewald  sphere  construction  yields  a relationship 


where  \ is  the  electron  wavelength  in  angstroms,  g the  magnitude  of 

diffraction  spot  from  the  origin,  respectively.  If  the  diffraction 

agreement  with  Figure  2 (the  vector  s points  in  negative  y-direction) 
and  the  geometrical  relationship  expressed  in  equation  (2.10). 

In  most  cases,  deviation  parameters  are  reported  as  dimension- 
less quantities  by  multiplying  with  the  appropriate  two-beam  extinction 
distance: 

wg  - sg5g  (2.37) 


been  retained  also  in  the  case  of 


sphere  construction  can 


many-beam  calculations. 

If  only  systematic  reflections  are 


a relationship  f< 


the  deviation  parameter 


way  of  Figure  41)  in  a plane  to  find 
systematic  beam  ng  in  terms  of 


systematic  row  are  sufficient  to  fix  th 
Evald  sphere  with  respect  to  the  reciprocal  lattice  (since  Che  sphere 
must  always  pass  through  the  origin,  000).  The  calculation  of  further 

due  to  the  lack  of  symmetry.  The  derivation  carried  out  in  Appendix  A2 
is  based  upon  an  idea  due  to  Bell.  The  equations  have  been  pro- 
grammed for  solution  by  computer  (Appendix  A3.1). 

interest  can  in  principle  be  measured  on  the  diffraction  pattern. 
However,  in  order  to  avoid  calculating  profiles  for  unrealistic 
orientations  due  to  Inaccurate  measurements,  only  two  s values  should 


Coordinate  systems  describing  dislocation 
position  in  the  microscope,  y is  beam  direction. 

a positive  edge  and  a left-handed  screw  disloca- 


where  b and  b are  rhe  edge  and  screw  components,  respectively, 
render  the  logarithmic  argument 

will  somehow  alter  the  dislocation  strain  £ 
infinite  medium  conditions.  Unfortunately, 

information  about  the  exact  foil  surface 

at  the  appropriate  foil  surface.  However,  this 

displacement  equations  have  been  developed  by  Bullough  and  Newman'J 
describing  a screw  dislocation  in  a plate  with  either  traction-free 


rid  as  calculated  from 
10  general  analytical 
addielon,  there  is  no 
a practical 
literature  have 


with  c and  d being  Che  £oil  thickness  and  Che  depth  of  Che  dislocation, 
respectively. 


ysten,  i.e.,  as  a profile  crossing 
n x-direction.  The  displacement  fu 


s recorded  in  the  microscope 


[fly]  is  the  transformation  m« 
n the  x,y,z  system  and  the  u,' 


x (matiix  of  direction  cosines} 


(2.44) 


Furthermore,  the  origin  of  the  microscope  system  will  be  taken  at  the 
upper  foil  surface  with  the  dislocation  at  depth  d,  resulting  in  the 


, the  transformation  equations  r« 


2.4.S  Contrast  functions  g. | — j 


partial  derivatives  of  the  displacement  functions, 
the  components  of  these  contrast  functions  in  u,  v. 


(2 


Including  surface  effects,  the  expressions  for  a screw  dislocation 


In  equations  (2.47)  through  (2.49),  the  following  three  parameters 


defined: 


This  completes  the  theoretical  formulation  of  strain  contrast. 
A numerical  program  to  calculate  intensity  profiles  across  a general 
mixed  dislocation  with  infinite-medium  displacements  is  given  in 


Appendix  A3.1 


CHAPTER 


EXPERIMENTAL  PROCEDURES 


3.1  Intensity  Scanning  Device 

The  principal  effort  fror.  an  experimental  point  of  view  has 

The  apparatus  basically  Includes  two  components,  a detection 

altered  to  deflect  the  beam  in  any  direction  desired.  However,  for 

motor-controlled,  variable  resistor.  The  coil  support  has  been  designed 
cooling  requirements. 

multiplier  in  typical  scanning  applications  were  quite  similar  from 
size  of  a channeltron  electron  multiplier,  the  former  mode  of  operation 


Figure  4. 


offered  considerable  advantages  in  practice  since  the  entire  detection 
device  could  be  counted  on  the  small  viewing  screen  of  the  Philips 

through  binoculars  during  a scan  and,  in  addition,  to  remove  the  device 

circuitry  of  the  channeltron  detection  device  are  shown  schematically 


f distance  at 


■t  through  the  deflection  colls,  calibrated  in 
oss  the  image,  while  the  Y-axis  records  Che 

samples  the  intensity  distribution  carried  by  the  electron  beam  in  one 
particular  direction  which  has  conveniently  been  chosen  Co  be  per- 
pendicular Co  the  dislocation  line.  Magnification  and  resolution  of  c 
system  depend  on  operating  conditions  and  can  be  adjusted  Co  the 
requirements  of  a particular  application  by  selecting  microscope 
magnification,  aperture  site  and  X-Y  recorder  settings.  For  disloca- 
tion scanning,  profiles  were  recorded  on  a scale  ranging  from  about 

resolution  (independent  of  the  microscope  resolution)  of  5 8, 
Increasing  to  about  35  X at  the  lowest  magnification  which  has  been 

Within  Che  practical  range  of  dislocation  scanning  (at 
microscope  place  magnifications  of  10000  X and  above  corresponding  Co 


scale.  In  Figure  6,  Che  signal  curren 
Intercepted  by  the  small  viewing  screei 


Che  Kelchley 
electron  currenc 

linearity 


should  be  compared  to  the  typical  peak  height  of  dislocation  profiles 
as  Illustrated  In  Figure  6,  and  Is  seen  to  be  too  large  to  be  neglected. 
Although  Figure  6 could  In  principle  serve  as  a calibration 

fylnges) , this  has  not  been  attempted  in  the  case  of  dislocation 

establish  the  tero  intensity  level  for  each  profile.  Since  it  was 
In  many  cases  required  to  counter  image  instabilities  by  speedy  data 
recording,  Che  base  line  of  the  X-Y  recorder  was  often  shifted  to 
allow  many  profiles  to  be  recorded  without  changing  graph  paper.  In 

No  attempt  has  therefore  been  made  to  match  absolute  intensity 
levels  of  the  profiles  presented  in  Chapter  4.  However,  except  In 
particular  cases  (see  Chapter  4) , relative  peak  intensities  should 

Intensity  was  always  sdjusted  to  approximately  the  same  Keithley 
reading.  This  provided  one  of  the  reasons  why  mainly  bright  field 


Although  attempts  to  determine  tl 


ie  channeltron 


output  was  linear  wit 
a counting  technique 
Figure  6).  This  would  leave  only 
have  Co  be  a function  of  signal  ci 

Kelthley  electrometer  as  possible 
The  channeltron  gain  is,  < 

Profiles  should  therefore  not  be  i 

Although  a very  slow  decrease  in 


pulses  per  unit  of  time  (< 
' the  channeltron  gain  (wh: 
turrenc  in  contradiction  tt 


-ever , time-dependent,  particularly 

corded  during  Che  first  fifteen 

in  was  still  observed  after  two 
effect  will  not  be  noticeable 


3.2  Sample  Preparation 

Most  of  Che  experimental  work  was  carried  out  on  aluminum  foils 
obtained  by  standard  evaporating  techniques.  The  aluminum  was  deposited 
at  approximately  350*C  and  10"°  Torr  on  100  surfaces  of  single-crystal 

After  evaporation,  a foil  was  cut  on  Che  substrate  into  small  square 
sections  which  were  then  floated  on  water  and  transferred  to  standard 
200  mesh  microscope  grids. 

<111)  in  the  foil  plane,  and  exhibited  irregular  dislocation  networks 


In  addition  to  the  vapor-deposited  films,  some  thin  foils 
were  prepared  by  electrolytic  chinning  of  bulk  aluminum.  They  were 
taken  from  a 1/4  inch  diameter  single-crystal  bar  of  high-purity  aluminum 
which  was  deformed  about  2 percent  in  torsion  Co  increase  the  probability 

diameter  discs  were  eleccrolytically  jet-dimpled  and  polished  until 
the  first  hole  appeared. 

The  same  procedure  was  employed  to  obtain  Cu-Al  foils  from  a 
grown  with  a 4 by  4 on  cross-section.  These  samples  were  used  to 


record  some  profiles  across  stacking  faults  and  extended  dislocations, 
some  of  which  are  presented  in  Section  4.S.  The  polishing  conditions 
and  electrolyte  solutions  are  shown  in  Table  2. 


3.3  Contrast  Parameter  Measurements 

evaluate  foil  orientation,  foil  thickness  and  dislocation  depth  for 
a particular  contrast  experiment,  where  emphasis  is  placed  on  the 


A diffraction  pattern  carries  all  the  information  necessary  t< 
of  the  Kikuchi  line  pattern  produced  by  the  inelastically  scattered 


L-defined  Klkuchl  lines 


if  the  electron  beam, 
can  be  caused  by  both  plascic  ar 

patcerns  In  Figure  7a  end  7b  is 


elastic  deformation  i 
The  difference  between 
due  to  the  larger  amou 


.1  bending  associated  with  a larger  diffracting  area, 
'e  become  very  hazy  in  Figure  7c  due  to  the  presence 


Klkuchl  lines 


It  is  apparent  from  equation  (: 
distance  between  the  diffraction  spot  • 

fraction  of  the  distance  of  the  diffraction  s; 
le  quantity  to  be  measured  on  the  pattern,  cal! 

! pattern  symmetry.  This  method 

surucy  depends  upon  the  reciprocal  intensity  spikes  being  perpendicular 
the  plane  of  the  diffraction  pattern.  For  this  reason,  most  of  the 
itrast  experiments  have  been  carried  out  close  to  the  111  orientation 


by  using 


The  exact  position  of  the  Klkuchl  line  of  Interest  can  be 
measured  more  accurately  from  higher-order  lines  which  are  always 
sharper,  since  s.  is  proportional  to  s /g  . Thus,  in  Figure  8,  the 


half-way  between  220  lines.  The  change  in  position  of  a 220  Klkuchl 
line  for  a variation  of  1.0  Id  V is  indicated  in  Figure  8 and  clearly 
illustrates  that  it  is  rather  difficult  to  attempt  to  measure  w within 


(C) 


Figure  7.  Influence  of  crystal  deformation  on  Kikuchl  patterns. 
Dislocation-free  grain  in  (a)  and  (b)  with  5li  and 
25y  diffraction  aperture,  respectively.  Adjacent 
grain  in  (c)  containing  many  dislocations,  5p 
aperture. 


Figure  8.  Determination  of  s.  from  diffraction  pattern.  A 
change  of  1.0  in  w causes  the  220  Kikuchi  line  to 
shift  by  flSj.  Beam  direction  close  to  111. 


S3 


orientation),  Klkuchl  lines  were  much 


.1-deflned  Klkuchl  pattern. 

t of  the  Klkuchl  line  for  a given 

im  the  plane  of  a foil  In  111 
•e  diffuse  than  those  generally 


because  it  yields  fails  of  constant  thickness  wl 
by  optical  interferometry,  allowing  for  a nominal  accuracy  of  50  X. 

In  practice,  foil  bending  on  the  substrate  produced  curved  Interference 
fringes,  and  thickness  values  were  probably  uncertain  within  about 
200  8,  corresponding  to  an  accuracy  of  i 32  for  a typical  3000  8 foil. 

were  In  agreement  with  the  number  of  grain  boundary  fringes  counted 
on  electron  micrographs  which,  however,  allow  the  foil  thickness  In 
aluminum  to  be  determined  within  only  500  8 at  best,  using  the  111 
reflection. 

by  measuring  the  width  of  slip  traces  purposely  Introduced  into  some 
e (see  Figure  9)  led  to  thickness  values  smaller  than  those 


methods. 


seen  moving  slowly  along  the  slip  traces  when  the  foil  was  nc 
strained  (dislocations  at  A,  B and  C in  Figure  9).  Taking  tl 

width,  good  agreement  was  found  with  interferometry  data. 

depth,  Intensity  scans  were  originally  planned  to  be  taken  at  the 
center  of  long  straight  dislocations  sloping  through  the  foil.  Howe\ 
in  both  vapor-deposited  and  electrolytically  thinned  specimens,  these 

for  a perpendicular  orientation  with  respect  to  the  foil  surface. 
Figure  10  shows  how  this  problem,  at  least  in  some  cases,  could  be 
intersection  of  dislocations  with 

e straight  dislocation  in  Figure  1( 

thickness.  A contrast  experiment  with  a 220  reflection 

At  this  point  of  the  discussion,  it  should  be  e\ 
is  a function  of  the  operating  reflection  causing  concrg 


s reference  for  i 


distance 

) X 

, - 1036  X) 


experiments  ro  be  presented  in  the  following  chapter,  220  reflections 

previously,  it  minimized  the  error  introduced  into  the  precise  deter- 
mination of  the  center  of  diffraction  spots  which  appeared  to  be  most 
critical  in  obtaining  accurate  values  for  the  deviation  parameter. 


DISCUSSION 


i-%  - 2 


directions  of  Burgers  vectors.  This  procedure  proved  to  be  reliable 
and,  for  foils  parallel  to  the  slip  plane,  much  faster  than  the  usual 
g*b  " 0 determination.  In  addition,  it  often  provided  1: 
information  about  the  dislocation  type  during  observation. 


applications  and  will  be  included  in  the  present  section  on  results. 

if  g'b  - 2 (two-beam  case).  This  fact  may  be  spplied  to  advantage  in 

Will  then  exist  only  one  220  reflection  which  can  give  g.b  - 2 for 
a particular  dislocation,  and  this  220  reflection  must  be  parallel  to 
the  Burgers  vector.  Thus,  a single  two-beam  orientation,  where  a 220 
diffraction  vector  gives  rise  to  a double  image,  can  uniquely  determine 
the  direction  of  b. 

Full  advantage  of  such  a procedure  can  be  enjoyed  when  Che 
dislocation  slip  plana  coincides  with  or  is  near  to  the  foil  plane. 


reflections  which,  at  zero  tilt,  are  already  very  close  to  the  Ewald 
sphere  (111  diffraction  pattern).  A small-angle  tilt  of  the  specimen 
must  therefore  lead  to  a double  image  with  g 9 220  and  thus  establish 


the  direction  of  b.  In  addition,  since  Che  dislocation  line  is  not  far 
removed  from  the  image  plane,  the  relative  rotation  between  diffraction 

edge  or  mixed).  Sample  rotation  need  only  be  applied  if  the  particular 
220  vector  is  very  close  to  the  tilt  axis,  or  if  the  220  reflection 
fails  to  produce  reasonable  two-beam  contrast  while  passing  through  the 
rarely  been  encountered  in 


ones,  over  a wide  range  of  tilt,  which  produce  strong  contours,  and  can 
therefore  be  recognized  quite  easily  on  the  viewing  screen. 

Using  the  procedure  described  above.  Burgers  vectors  of 
dislocations  of  interest  could  be  found  within  minutes.  In  addition, 
a dislocation  of  known  character  could  easily  be  selected  for  detailed 
study.  For  this  purpose,  a 220  dark  field  contour  was  moved  over  a 
general  area  containing  dislocations,  and  the  line  direction  of  dis- 
locations having  double  images  was  compared  to  the  direction  of  the 
diffracting  vector. 


In  order  to  examine  the  reliability  of  the  g*b  • 2 criterion, 
numerous  b-determinations  were  carried  out  applying  both  the  g-b  - 2 
and  g.b  * 0 criteria,  and  agreement  was  found  invariably.  Furthermore, 
the  occurrence  of  a double  image  for  g.b  * 2 was  investigated  for  a 
variety  of  contrast  conditions  in  both  experiment  and  theory.  Figure 
11  shows  several  forms  of  appearance  which  double  images  may  exhibit 
in  a given  situation.  In  Figure  11a,  the  g*b  - 2 contour  crosses  a 
screw  dislocation,  and  it  is  evident  that  the  double  image  exists  only 


at  s-values  close  to  zero.  Figures  lib  and  11c  show  Inclined  mixed 
dislocations  sloping  through  Che  foil.  A long  edge  dislocation  in  the 
plane  of  the  foil  is  shown  in  Figure  lid  where  the  "doubling"  of  the 
double  image  is  typical  for  calculated  profiles  at  some  particular 
values  of  foil  thickness  (e.g.,  see  Figure  12  at  t ■ 2.8  £ ).  Figure 

another  potential  source  for  double  images.  Screw  dislocations  exhibit 
three-beam  double  images  (g.b  ■ 2 and  3)  at  A and  B,  and  two-beam 


Two-beam  profiles  with  g.S  ■ 2 and  v ■ 0 have  been  calculated 
for  many  different  values  of  foil  thickness  and  dislocation  depch. 

The  double  images  are  strictly  symmetrical  only  in  dark  field,  when 
the  dislocation  lies  at  the  center  of  the  foil.  In  this  case,  the  nature 
of  the  intensity  peaks  (bright  or  dark)  as  well  as  the  magnitude  of 

Figure  12.  The  change  of  image  contrast  with  thickness  in  Figure  12c 
is  very  similar  to  chat  predicted  by  the  calculated  profiles  in  Figures 

bottom  of  Figure  12c  have  been  used  to  Insure  both  v = 0 and  a knowledge 
of  t along  the  dislocation  line.  Although  the  image  at  t ■ i ( seems 
to  agree  better  with  an  absorption  parameter  of  B • 0.05  (Figure  12b) 
chan  with  B >0.10  (Figure  12a),  no  conclusion  can  be  drawn  since  the 

The  main  influence  of  varying  dislocation  depth  for  the 
g.b  ■ 2 profile  appears  as  a change  in  relative  heights  or  widths  of 


at  certain  combinations  of  dislocation  depth  and  foil  thickness  (see 
Figure  lib)  which  also  seem  to  be  affected  by  absorption.  The  g*b  * 2 
contrast  can  therefore  only  be  applied  as  an  affirmative  criterion  in 
the  Burgers  vector  determination,  i.e.,  it  cannot  be  used  in  order  to 
exclude  a particular  220  reflection  from  giving  the  direction  of  b 
before  the  double  image  has  been  found. 

Finally,  as  with  all  criteria  for  the  determination  of  Burgers 

double  images  of  edge  and  screw  dislocations  have  also  been  reported 

which  are  known  to  be  strong  in  this  material.  On  the  other  hand, 

tion  which  is  steeply  inclined  with  respect  to  the  diffraction  plane, 

Although  none  of  these  effects  would  apply  specifically  to  the  contrast 
produced  by  a 220  reflection  for  long  dislocations  in  a 111  foil 
orientation,  contrasc  anomalies  in  particular  cases  can  never  be 

in  the  history  of  the  g.b  ■ 0 criterion. 


t successful  applications  o: 


the  dynamical  t' 


I.w  - 1. 


1 surface.*6*  The  exact  nature  of  the  oscillations 
was  found  to  depend  on  bright  field  or  dark  field  conditions  as  well 
as  on  the  foil  thickness,  and  the  pertinent  extinction  distance 
appeared  to  give  the  period  of  oscillations.*6'39*  Figure  13  shows 
the  basic  features  of  these  periodic  profile  changes  with  dislocation 

essentially  identical  profile  shape  is  obtained  when  the  dislocation 
depth  position  has  changed  by  one  effective  extinction  distance,  which 
s simply  given  by  S (Figures  13a  and  13b),  but  which  changes 
for  w i 0 (Figure  13c) . Deviations  from  this 
periodicity,  however,  may  occur  when  the  dislocation  approaches  the  foil 

The  profiles  in  Figure  13  are  calculated  for  depth  positions  at 
every  quarter  of  the  repeat  distance.  Additional  computations  have 
shown  that  the  profile  change  is  gradual  in  between,  in  good  qualitative 
agreement  with  the  experimental  images  of  Inclined  dislocations  in 

Almost  identical  profiles  to  those  presented  in  Figure  13  have 
been  found  for  g*b  * 1 when  varying  the  foil  thickness  while  keeping 

this  case  showed  a period  of  2£  for  the  foil  thickness,  again  corres- 
ponding to  a change  of  one  f,  in  dislocation  depth.  It  is  therefore 
expected  that  variations  in  either  thickness  and  depth  position  lead 
to  equivalent  changes  in  the  general  shape  of  the  dislocation  profile, 
depending  only  on  the  relative  positions  of  dislocation  and  exit  surface 


Oscillatory  contrast  in  e dark  field  Image 

(g-S  ■ 1 and  v = 0)  of  mined  dislocations  inclined 


with  respect  to  the  periodic  intensity  distribution  in  the  foil  which 
gives  rise  to  the  familiar  thickness  extinction  fringes. 

For  the  purpose  of  image  profile  matching,  it  may  be  concluded 
that  profile  changes  are  to  be  expected  over  a range  of  thickness  and 

a contrast  experiment.  The  same  situation  has  been  found  for  small 
variations  in  the  deviation  parameter  w where,  however,  profile  changes 

rules  for  profile  matching  can  therefore  be  established: 

(1)  In  order  to  distinguish  between  two  possible  choices  for 
the  correct  displacement  field,  profiles  should  be  cal- 


Foil  thickness  or  dislocation  depth  need  i 


For  a definite 


A large  number  of  calculations  for  g.S  = 


revealed  a quite 


Co  determine  Bj^j  in  e copper-silicon  alloy.  In  Che  following,  present 
resales  for  Che  111  reflection  in  aluminum  will  be  described. 

Ic  has  been  found  that  g*o  ■ 1 images  in  either  bright  field  or 

effect  on  the  profiles  of  type  El  and  ZZZ  (only  relative  Intensity  is 

way  as  an  Increase  in  foil  thickness,  will  move  Che  image  peak  in 
profile  type  ZZ  closer  to  the  true  position  of  the  dislocation  line, 
and  will  gradually  suppress  the  bright  intensity  peak  in  profile  type 

profile  type  ZZZ  represents  a particularly  suitable  application  for 
profile  matching.  Unfortunately,  representative  images  could  only  be 

method.  However,  as  shown  below,  qualitative  image  interpretation 
tal  the  existence  of  an  apparent  anomalous 
:h  decreases  with  increasing  thickness. 

:s  some  selected  black-white  images  of  inclined 
iminum,  showing  clearly  i decrease  in  intensity 
'oil  thickness.  This  effect  agrees  with 


absorption  parameter  wl 


calculated  profiles  of  type  111,  some  of  which  are  given  in  Figure  16, 


Black-white  images  (profile  types  1 end  III)  of 
inclined  mixed  dislocations  in  aluminum  at  w » 0 
as  a function  of  foil  thickness.  Dark  field  in 
(a),  (d)  and  (e),  bright  field  in  (b)  and  (c). 
Two-beam  diffraction  vector  g B 111,  scale  500  S. 


Cion  for  g-b  = 1 and  v = 0 as  a function  of  the  absorpeion 

are,  respeccively:  3.0/1.75  (a),  5.5/2.75  (b) , 8.0/4.25  (c), 
and  11.0/5.75  (d).  Dark  field  in  (a),  (c)  and  (d),  bright 
field  in  Cb) . Intensity  level  is  indicated  for  0.01. 


different  thickness  values  but  Identical  absorption  parameters 
(diagonal  rows).  The  profiles  have  been  calculated  for  a screw 
dislocation  after  only  a very  small  difference  In  peak  width  was  found 
In  profiles  computed  for  the  mixed  dislocations  In  figure  15. 

Closer  inspection  will  show  that  the  intensity  of  white  lobes 
in  Figure  15  does  not  decrease  as  rapidly  as  predicted  by  calculations 
involving  constant  B.  For  example,  a minimum  absorption  parameter  of 
B ■ 0.07  in  Figure  16a  is  necessary  to  produce  visible  dark  lobes  at 

Figure  16c  leads  to  a very  weak  bright  peak  at  t - 8 5 which  does  not 
correspond  to  the  rather  bright  lobes  in  Figure  15d.  At  t - 11.0  C , 

B • 0.07  would  all  but  suppress  the  bright  peak  completely,  whereas  the 
still  visible  white  lobes  in  Figure  15e  (compare  black  and  white  lobes 
along  the  dislocation  line)  show  better  agreement  with  an  absorption 


Although  this  qualitative  comparison  between  images  and  profiles 
can  only  give  a general  indication  that  the  apparent  ai 
absorption  parameter  B tends  to  decrease  with  increasing  f< 

from  thickness  fringes.  Applying  equation  (2.34)  together  with  equation 
(2.35)  to  an  evaluation  of  intensity  traces  across  thickness  fringes 
in  aluminum  and  copper-aluminum,  a consistent  initial  decrease  in  B 
been  observed  in  thin  portions  of  Che  foils. 


example  for  this  behavior  will  be  given  below  (Figure  18) . In  addition, 
a difference  in  Che  measured  value  of  B has  been  observed  by 


Hashimoto 


,(26) 


,(42) 


ie  measured 


thick  foils,  and  Takagi  and  Ii 
te  correct  anomalous  absorption  parameter  mi 

face  already  suggested  such  an  approach  after  an  investigation  of 

these  additional  electrons  in  the  imaging  beam  is  well  documented  by 
the  dependence  of  measured  absorption  parameters  on  the  sire  of  the 
objective  aperture^®* 27,43)  an^  on  tjie  f0ii  temperature.^**^ 

In  order  to  test  the  dependence  of  the  apparent  anomalous 

linearity  with  the  help  of  the  calibration  curve  in  Figure  6,  the  B-values 
plotted  as  function  of  foil  thickness  in  Figure  18  have  been  calculated 
from  equations  (2.34)  and  (2.35).  The  decrease  of  B with  increasing 


thickness,  as  discussed  above.  Is  evident.  On  the  other  hand,  the 
circles  with  arrows  in  Figure  18  roughly  Indicate  the  range  for  B as 
estimated  from  comparison  between  the  black-white  images  of  Figure  15 
and  the  calculated  Intensity  profiles  in  Figure  16.  It  should  be 
noticed  that  the  dislocation  in  Figure  15d  (as  well  as  the  dislocation 
in  Figure  15e)  is  slightly  displaced  from  the  center  of  the  dark  thick- 
ness fringe,  thus  somewhat  reducing  the  intensity  of  the  bright  lobes 
(there  are  no  white  lobes  when  the  dislocation  lies  between  dark 
fringes,  see  Figure  19),  and  probably  leading  to  an  overestimate  of  B 


snee  from  dislocation  contrast  analysis 
satisfactory,  and  an  apparent  ai 


absorption, 

absorption  parameter  as  obtained  from  thickness  fringes  appears  t( 
promise  a successful  natch  between  theoretical  and  experimental 
dislocation  profiles.  In  addition,  comparison  of  the  solid  and  bi 
curves  in  Figure  18  with  the  encircled  estimates  of  an  empirical  I 


choice,  and  it  Indeed  shows  better  agreement  with  theoretical  calculations 
of  B at  large  foil  thickness  (see  Table  1 in  Section  2.3). 


Dislocation  1: 


The  oscillating  images  of  dislocations  in  Figure  19,  repre- 
senting profile  types  1 and  11,  clearly  show  chat  Che  position  of  the 
dislocation  image,  even  at  w - 0,  does  not  in  general  coincide  with 

thick  foils  (Figure  19e).  Although  difficult  to  measure  accurately, 

increasing  foil  thickness  in  Figure  19  is  somewhat  slowed  by  a simul- 
taneous decrease  in  the  apparent  anomalous  absorption  parameter. 

In  order  to  test  whether  the  dynamical  theory  predicts  a 

le  dislocation  image  with  respect  to  the  disloca- 
te the  mixed  dislocation  in  Figure  19d  at  t ■ 8 £ 


profiles  by  cutting  off  their  peaks  to  obtain  tl 
amplitude  of  the  oscillations  with  respect  to  tl 


■ 0.05.  In  Figure 


dislocation  line, 


Oscillatory  ineges  (profile  types  1 and  II)  of 

as  a function  of  foil  thickness.  Dark  field  in 
(a),  (c)  and  (c),  bright  field  in  (b)  and  (d). 


Bright  field  of  inclined  mixed  dislocation  in 
aluminum  at  g-S  - 1 and  w • 0 (profile  type  1 and 
II).  Foil  thickness  t - 8 £ . Braving  indicates 
oscillations  predicted  by  calculations.  See  text 


Figure  21  shows  a grain  containing  long  dislocations  in  almost 
perfect  screw  orientation  which,  of  many  dislocations  studied  in  the 

!St  suited  for  a matching  experiment.  First,  surface 
iken  into  account  for  a screw  dislocation.  Secondly, 


spacing  of  the  optical  interferometer  whi 
accuracy  of  measuring  t.  Finally,  it  cot 
grain  boundary  fringes  (see  also  Figure  3 

location,  a 220  contour  was  chosen  which 
dislocation  line.  In  additicn,  this  reflection  gave  g.b  - 2,  and  it 
had  been  observed  previously  that  g.S  - 2 produces  stronger  contrast  i 
large  values  of  w than  does  g.S  » 1.  The  intensity  profiles  measured 
at  point  C (Figure  21)  for  six  different  orientations  are  shown  in 
Figures  22b  and  23b,  where  the  « 


to  be  equal  to  the  reference 

3 be  established  from  the 
) that  the  dislocation  marked 
vicinity  of  point  C. 

in  nearly  parallel  to  the 


from  the  diffraction  patterns  in  Figure  21.  The  theoretical  profiles 
calculated  for  the  experimental  set  of  contrast  parameters  are  given  in 
Figures  22a  and  23a  (g-b  ■ +2  for  all  six  profiles),  and  Figures  22c 
and  23c  show  computed  profiles  which  gave  a better  match  with  experiment. 

Most  characteristic  for  the  dislocation  images  represented  in 
Figures  22  and  23  appears  to  be  the  intensity  change  from  main  peak  to 
background,  and  it  can  be  seen  that  good  agreement  is  obtained  when  an 


Screw  dislocations  in  bright  field.  The  arrow  at 
C narks  position  where  the  scans  of  figures  22  and 
23  were  token  at  the  foil  orientations  represented 
by  the  diffraction  patterns  above. 


per  division.  Foil  thickness  t “ 2.8  Cg.  See  text  for 


accuracy  limit  of  ± 0.1  for  w is  assumed.  Such  error  limit  is  actually 

diffraction  patterns  in  Figure  21.  Thus,  even  better  agreement  between 
theory  and  experiment  may  be  obtained  simply  by  calculating  additional 
profiles  for  a larger  and  more  extensive  range  of  w-values.  In  this 
connection,  it  should  be  emphasized  that,  at  present,  profile  matching 
in  diffraction  contrast  is  restricted  to  a trial-and-error  procedure  in 
which  the  degree  of  agreement  depends  somewhat  on  the  number  of  computed 


4.3.2  Sense  of  Burgers  vector 

It  can  also  be  seen  from  Figures  21,  22,  and  23  that  the  calculated 
curves  immediately  match  the  experimental  profiles  taken  with  the  022 
reflection  (w  ■ 0.8  and  w ■ 1.7),  whereas  they  must  be  inverted 
(equivalent  to  changing  g.6  - +1  to  g.S  - -1)  in  order  to  match  profiles 
produced  with  the  022  reflection.  The  positive  Burgers  vector  in 
Figure  21  is  thus  to  the  left  (opposite  to  the  positive  line  direction 

4.3.3  Relative  intensities 


:ar  nothing  hi 


e scanning  system  wi 


due  to  relative  Image  rotations  at  different  foil  orientations  (see 
diffraction  patterns  in  Figure  21),  the  horizontal  scale  was  not  the 
same  for  all  scans  (the  x-scale  depends  on  the  scanning  direction,  see 
Section  3.1).  The  profiles  in  Figures  22  and  23  have  therefore  been 


enlarged  photographically  by  a factor  ranging  from  1.15  to  1.35  in 
order  to  compare  measured  peak  width  to  calculations,  and  the  original 
peak  heights,  in  terms  of  the  respective  calculated  peak  intensities, 


k height  calculated 


■tter  in  relative  peak  Intensities  si 
e nonlinearity  of  the  scanning  systi 
the  fact  that  the  diffraction  pattei 


A more  probable 


beam.  Thus,  the  profiles  i 
A. 3. A Surface  effects 


Figures  22  and  23  may  have  been  measured 
s with  respect  to  calculations.  This 


dislocations,  surface  effects  on  the  image  of  the  screw  dislocation  in 
Figure  21  should  be  investigated,  Matching  in  Figures  22  and  23  was 
achieved  with  profiles  calculated  from  the  displacement  equations  (2.48) 

with  rigid  surface  conditions  have  been  employed  to  calculate  profiles 

w for  the  two  experimental  scans  in  the  figure.  The  main  difference 
from  Che  case  of  free  surfaces  appears  to  be  the  persistency  of  the 
double  image  for  g,b  • 2 for  larger  values  of  the  deviation  parameter 
under  rigid  surface  conditions.  This  is  qualitatively  in  contrast  to 
observation  on  the  microscope  screen  where  the  double  images  disappeared 
immediately  when  tilting  away  from  w - 0 (see  also  Figure  11a).  From 
a matching  point  of  view,  free  surface  conditions  (Figures  22  and  23) 
clearly  provide  better  agreement  than  the  presence  of  a rigid  surface. 

in  aluminum  disappears  quickly  due  to  stress  relaxation  at  the  surface/9' 
since  stress  relaxation  would  not  occur  easily  in  the  presence  of  an 
elastically  rigid  surface  layer. 

The  experimental  evidence  favoring  the  free  surface  conditions 
for  dislocation  displacements  is  quite  welcome  for  the  calculation  of 
intensity  profiles.  It  has  been  shown  previously*47'  that,  as  long  as 
the  dislocation  lies  at  a distance  of  5 from  the  foil  surface,  its  main 
image  peak,  as  calculated  from  infinite  medium  displacements,  is  not 
affected  by  taking  a traction-free  surface  into  account,  the  only 
difference  being  a faster  approach  of  the  profile  to  background 
intensity  at  some  distance  from  the  main  peak  due  to  stress  relaxation. 


Consequently,  Infinite  medium  displacements  could  be  employed  to 
calculate  Intensity  profiles  of  edge  or  mixed  dislocations,  without 
altering  the  matching  conditions. 


Hith  the  exception  of  the  profile  for  w ■ 1.7,  all  experimental 
images  in  Figures  22  and  23  are  broader  than  predicted  by  calculations. 
Peak  broadening  may  be  expected  to  occur  due  to  the  divergence  in  the 
incident  beam,  or  due  to  diffuse  intensity  produced  by  inelastic 
scattering.  In  both  cases,  the  degree  of  broadening  should  depend 
upon  the  choice  of  the  operating  reflection  (see  Section  4.3.6  for 
detailed  discussion).  Ill  and  311  profiles  have  therefore  been  matched 
to  study  the  relative  peak  width  as  a function  of  g. 

Figure  25a  shows  a mixed  dislocation  imaged  in  bright  field 
with  a 311  reflection,  from  which  the  experimental  scans  in  Figure  26 
have  been  taken.  The  curved  dislocation  in  Figure  25b  has  supplied  the 
111  profiles  shown  in  Figure  27.  In  both  examples,  information  about 

For  the  dislocation  in  Figure  25a,  the  depth  could  not  be  determined. 
Profiles  were  therefore  calculated  for  various  depth  positions,  and 
it  was  found  that  matching  within  the  accuracy  of  Figure  26  could  be 
achieved  for  more  than  one  depth  position,  provided  the  deviation 
parameters  were  adjusted  accordingly.  If,  however,  the  dislocation 
was  assumed  to  lie  at  or  near  the  center  of  the  foil  (calculated 
profiles  in  Figure  26),  the  match  was  obtained  over  a range  of  w-valuea 
which  are  in  excellent  agreement  with  the  diffraction  patterns  in 


From  boundary  fringes  of  the  grain 

the  left  end  of  the  dislocation  was  taken  ii 
period  of  the  effective  extinction  distance, 
parameter  in  this  case  col 


the  deviation 

s case  could  only  be  estimated  roughly  as  w < 1.0,  due 
ill-defined  Klkuchi  lines.  From  profiles  calculated 
taking  into  account  Che  different  nature  of  the 
le  scanning  positions  due  to 


5te  that  the  w-values 

Lso  note  the  intensity 

calculated  profiles  is 

distance, 
different 


dislocation  at  the 

match  was  obtained  as  shown  in  Figure  27. 
reflect  the  experimental  observation  that 
away  from  the  dislocation  at  C than  at  A. 
scale  which  is  magnified  10X  with  respect 

due  to  a larger  foil  thickness  in  terms  of 
Such  a large  intensity  difference  between 
reflections  has  not  been  observed  experimentally.  This  suggests  th 
mean  absorption  parameter,  A,  may  be  less  chan  the  value  of  0.1  whl 
has  been  assumed  in  calculations. 
d.3.6  Relative  peak  width 

In  Figure  28,  relative  peak  widths  defined  as  the  ratios  be 
measured  and  calculated  half  widths  have  been  plotted  as  a function 
(h  + k + 1 ) for  the  three  reflections  whose  dislocation  images  h. 
been  matched  in  the  previous  sections.  The  data  can  only  be  taken  j 
a rough  indication  for  the  peak  widths  because  it  is  difficult  to 
h peak  heights  of  the  experimental  profiles  objectively. 


Figures  22,  23,  26  and  27,  as  a function  of 
diffraction  vector  g. 


The  large  scatter  In  profile  width  for  the  220  and  311  reflec- 
tions nay  be  due  to  Imperfect  Image  matching,  but  could  also  be  explained 
by  an  image  drift  in  the  microscope  which  was  often  observed  during 
contrast  experiments  (scans  were  discounted  when  this  was  noticed). 

appreciably  affected  by  the  nonlinear  Intensity  scale  of  Che  scanning 
system  because  all  profiles  were  taken  at  approximately  the  same  level 
of  output  intensity  (see  Figure  6). 

Despite  the  possible  errors  involved  in  the  data  of  Figure  28, 

1)  the  experimental  profiles  tend  to  be  broader  than  predicted 

2)  the  degree  of  broadening  seems  to  increase  with  the 
magnitude  of  the  diffraction  vector, 

4.3.7  Image  broadening 

Incident  beam  divergence  constitutes  a possible  source  for 
image  broadening  because  It  will  lead  to  an  effective  range  of  w-values 
over  which  a given  set  of  planes  will  diffract  electrons  In  the  foil. 

If  Che  Incident  beam  deviates  by  48  from  the  Sragg  direction,  then 


"-«g|g|48  (4.1) 

diffraction  takes  place  at  a given  beam  divergence  will  Increase  with 
1*1"  “here  n must  be  larger  than  one.  Although  the  effect  of  a Aw  on 


the  image  width  will  probably  be  a function  of  ti 
conditions  (i.e.,  of  the  degree  of  variation  witl 
profile  exhibits),  broadening  should  in  general 
for  a larger  <lw.  Beam  divergence,  which  has  necessarily  fc 
rather  large  in  order  to  provide  for  sufficient  intensity 
s explain  the  data  in  Figure  2£ 


particular  diffraction 


h fully  focussed  ai 


detection  system,  could  thu 

the  same  dislocation  image  produced  under  b< 
defocussed  illumination.  Figure  29  compares  the  results.  A 60  sec 
exposure  of  the  weak  image  in  Figure  29a  still  produced  less  blackening 
on  the  plate  than  a 4 sec  exposure  for  the  high-intensity  image  in 


n both  cases,  the  dislocation  hi 
e Intensity  distribution  across 


verified  by  tl 

effect  of  beam  divergence  on  the  peak  width  sho 

tion  for  positive  and  negative  deviations  in  w. 

Despite  these  rather  extreme  conditions 
is  noticed  in  Che  profiles  shown  in  Figure  29. 

fact  chat  the  plate  was  slightly  underexposed. 


In  addition,  the 
rue  dislocation  posi- 

ic  apparently  wider 
ie  explained  by  the 
ils  would  essentially 


be  undertaken  before  final  judgement 


The  second  potential  source  for  image  broadening,  the  diffuse 
intensity  due  to  inelastic  scattering,  is  also  expected  to  be  more 
effective  for  larger  diffraction  vectors:  The  dynamically  scattered 
part  of  the  electron  wave,  i.e. , the  diffracted  intensity,  will  decrease 
with  increasing  g,  while  the  diffuse  intensity,  to  a first  approximation, 
remains  constant.  Since  the  Bragg  spot  on  a diffraction  pattern  is 
surrounded  by  a rather  large  area  of  diffuse  intensity,  the  image  width 
should  in  this  case  depend  upon  Che  site  of  the  objective  aperture. 

Figure  30  presents  the  results  of  a contrast  experiment  where 
both  bright  field  (Figure  30a)  and  dark  field  (Figure  30b)  profiles 


have  been  taken  with  objective  apertures  of  5u 
irregularities  produced  by  Intensity  noise,  the 
heights  due  to  the  larger  aperture  in  bright  fi 
noted.  This  reduction  is  not  observed  (within 
fluctuations)  in  dark  field,  and  can  therefore 
effect  of  lnelastlcally  scattered  electrons,  th 
centered  around  the  forward  beam  direction.  Th 

Intensity  to  fill  up  the  intensity  valleys,  and  coi 
peak  broadening. 

Finally,  a comment  should  be  made  about  thi 
electrons  which  have  lose  only  a small  fraction  of 
inelastic  collisions  in  the  foil.  It  has  been  shot 
trons,  to  some  degree,  follow  diffraction  laws  and 
tion  about  a defect.  <4®_54)  They  are  thought  to  hat 
strongly  to  the  dark  field  profile  in  Figure  31  whi 


Id  (Figure  30a)  should  1 
lie  limit  of  intensity 
e Interpreted  as  an 
majority  of  which  are 
reduction  in  peak 
positive  peak 


: contributed 
h was  originally 


Experimental  profiles  across  a 
w = 0,  under  bright  field  cond 


expected  to  show  image  broadening  to  a lesser  degree  chan  e 
Che  case  of  bright  field  profiles.  However,  Che  broadening  effect 
in  Figure  31  is  actually  quite  pronounced.  On  Che  ocher  hand,  the 
experimental  profile  in  Figure  31  was  found  Co  exhibit,  in  comparison 
to  profiles  of  the  same  dislocation  with  w close  Co  aero,  more  than 
twice  the  intensity  predicted  by  calculation.  Such  an  increase  may  well 
be  explained  by  the  large  portion  of  inelastically  scattered  intensity 
which  oust  have  passed  through  the  objective  aperture  (see  Che  diffrac- 
tion pattern  in  Figure  31),  and  may  be  responsible  for  the  large  degree 

It  must  be  added  that  the  calculated  profile  in  Figure  31  was 
quite  Insensitive  to  the  dislocation  position  in  the  foil,  so  chat  Che 
lack  of  information  about  the  dislocation  depth  did  not  interfere  with 
. the  matching  procedure. 

4*3.8  Influence  of  dislocation  type 

If  matching  is  to  be  applied  successfully  to  identify  unknown 
displacement  fields,  it  must  be  possible  to  distinguish  screw  or  edge 
character  of  a given  dislocation  profile.  As  long  as  the  slip  plane  is 
not  too  Steeply  inclined  with  respect  to  the  diffraction  plane  (<«•), 
the  edge  displacement  component  perpendicular  to  the  slip  plane  will 
have  little  effect  on  the  image,  and  the  dislocation  character  is 
predicted  to  affect  only  the  profile  width. <6>  If  imaged  under 
identical  contrast  conditions,  the  edge  dislocations  should  exhibit  a 
profile  twice  as  wide  as  the  screw  dislocation.  The  experimental 
profiles  in  Figure  32  are  thus  in  reasonable  agreement  with  theory, 


Experimental  dark  field  profiles  approximate 
screw  and  edge  dislocations  for  g.b  = 1 (a)  and 
g*b  ■ 2 (b).  Scan  positions  indicated  by  arrows. 


for  both  g-b  - 1 (Figure  32a)  and  g-b  - 2 (Figure  32b) . Image 
contrast  in  Figure  32a  shows  a strong  Influence  of  dislocation  depth 
which  Is  typical  for  relatively  thin  foils,  and  the  similarity  of  the 
contrast  oscillation  at  the  two  scanning  positions  suggested  that 
profiles  here  would  be  taken  at  equivalent  depth  positions. 

When  deducing  the  dislocation  character  from  the  profile  width, 
it  must  be  remembered  that  mixed  dislocations  give  rise  to  profiles 
which  may  be  narrower  than  the  screw  image,  but  which  may  also  be 
wider  than  the  edge  profile/6’55*  Therefore,  a given  profile  width 
cannot  uniquely  determine  the  nature  of  the  dislocation.  It  is, 
furthermore,  obvious  that  image  broadening  as  discussed  above  is  an 
important  factor  in  such  a determination  of  displacement  fields. 

4.4  Manv-Beam  Effects 
4.4.1  Systematic  reflections 


Of  particular  interest  to  image  analysis  are  the  effects  of  the 
so-called  systematic  reflections  which  are  always  excited  to  some 
extent.  In  addition,  their  deviation  parameters  are  well  known  for  any 
particular  two-beam  orientation  end  can  be  calculated  from  equation 
(2.38).  A four-beam  calculation  for  the  111  reciprocal  lattice  row 
reported  by  Howie  and  Basinskl(16)  in  an  aluminum  foil  of  t - 8 { did 
not  show  significant  changes  with  respect  to  the  two-beam  profile. 


numerical  integration  of  equations  (2.30)  could  no  longer  be  carried 
out  within  the  limit  of  the  fifth-order  error  evaluation  employed  in 


became  too  large.  A similar  effect  was  reported  by  Thomas  and  Bell 
who  interpreted  equations  (2.30)  as  being  of  "nonconservacive"  nature. 
Although  the  difficulty  reportedly  can  be  avoided by  using  equation 

(2.27)  directly  for  intensity  calculations,  without  the  substitution 

(2.28)  to  remove  the  exponential  expressions  in  (2.27),  the  integrating 

After  removing  an  error  evaluation  altogether  in  the  Runge- 
Kutta  process,  systematic  four-beam  profiles  could  be  calculated  whose 
shape,  as  expected,  depended  upon  the  choice  of  step  size  for  the 
numerical  integration.  However,  such  a procedure  is  legitimate  if  the 
step  size  can  be  chosen  small  enough  so  thac  a further  decrease  will 
no  longer  alter  the  profile.  Figure  33  shows  some  of  the  results 
obtained  for  the  111  systematic  row  in  a thin  aluminum  foil  of  t - 3 £ , 
where  profiles  for  both  bright  field  and  dark  field  conditions  are  shown 
as  a function  of  integration  step  size,  Ay.  As  is  expected  from 
, the  influence  of  step  size  is  a function  of  the 
s well  as  of  Che  magnitude  of  s.  This  may  be  noted  in 
the  profiles  remain  essentially  constant  at  large 
le  dislocation,  but  depend  on  Ay  in  the  vicinity  of  the 

tep  size  of  0.01  £ , whereas  the  bright 
iws  a change  when  decreasing  Ay  to  0.005  £ . A similar 
ie  required  seep  size  was  found  in  the  choice  of  absorption 
Ay  ■ 0.02  £ , the  obviously  artifically  large  peak  in 
bright  field  as  well  as  Che  double  peak  in  dark  field  was  no  longer 
present  when  increasing  the  effect  of  absorption  by  changing  from 


displacement  R ai 
distances  from  cl 


correctly 


the  dislocation  line  result  in  erroneous  intensity  predictions  even  for 
the  two-beam  profile  at  w - 0 when  using  Ay  • 0.02  E . 

Without  the  integration  error  evaluation,  it  seems  therefore 
necessary  to  calculate  each  particular  image  profile  with  several 
Ay-values  in  order  to  exclude  artificial  features  introduced  by  the 
1 procedure.  Due  to  the  considerable  amount  of  computer  time 
e profiles  (Ay  - 0.005  Eg  required  more  than  3 minutes 
for  t - 3 Eg,  and  the  time  increases  almost  linearly  with  foil  thickness), 
the  study  of  systematic  many-beam  effects  did  not  advance  beyond  the 

Comparison  between  the  profiles  in  Figure  33  suggests  a strong 
e first  two  systematic  reflections  222  and  111  on  the 
:ase  with  w ■ 0,  in  this  case  mainly  represented  in  the  bright 
:e.  It  is  Interesting  to  note  the  apparent  disagreement  of 
•t  with  the  calculations  of  Howie  and  Basinskl^1^  mentioned 
.e  discrepancy  can  be  explained  by  the  much  thinner  foil  used 
e profiles  in  Figure  33  (t  ■ 3 c versus  t = 8 E ),  since 
it  is  known  that  the  intensity  on  the  reciprocal  relrods  at  a given 
distance  from  the  exact  diffraction  spot  in  general  decreases  with  foil 
thickness.  Thus,  higher-order  systematic  beams  are  expected  to  carry 
less  Intensity  in  thicker  foils,  an  effect  which  can  in  many  cases  be 
observed  in  the  microscope. 

The  experimental  two-beam  images  presented  in  the  preceding  section 
must  of  course  be  seen  in  the  light  of  these  first  results  on  systematic 
111  profiles  for  a foil  thickness  of 


involved  i: 


t - 7.90  Eg  should  be  those  predicted  by  two-been  theory,  other  profiles 
token  with  220  or  311  reflections  at  a foil  thickness  of  approximately 
3 C may  have  possibly  been  affected  by  the  presence  of  systematic  beams. 
On  the  other  hand,  the  larger  distance  between  spots  in  the  higher- 
order  reciprocal  lattice  rows,  combined  with  the  curvature  of  the  Ewald 
sphere,  leads  to  much  larger  deviation  parameters  for  the  systematic 
reflections  in  the  Utter  case  (a  factor  of  almost  3 between  s-values  for 


',  respectively).  In  fact,  the  four-beam  deviation 
was  apparently  large  enough  to  cause  totally 
s to  be  calculated  with. an  integration  step  size  of 


m effects  in  Figures  2\ 


1 reflection  in 


4.4.2  Konsvstematic  reflections 

It  should  be  noted  Initially  that  nonsystematic  reflections 
which  describe  all  diffraction  vectors  with  a direction  dl 
that  of  the  primary  two-beam  reflection  do  not,  in  principle,  : 

g.b  appears  in  the  Intensity  equations,  image  contrast  will  not. depend 
on  the  direction  of  g as  long  as  the  product  g*5  remains  constant. 

particularly  worthwhile  to  study;  first  the  magnitude  of  a third  beam 
effect  on  the  image  as  a function  of  s,  and  secondly  Che  feasibility  of 
matching  many-beam  Images. 

It  would  be  quite  desirable  to  know  how  far  a two-beam  orientation 


d (besides 


oust  be  improved  in  order  ro  really  represent  what  it  is  ca 
the  effects  of  systematic  reflections  of  course).  However, 
infinite  variety  of  possible  orientations,  it  is  probably  » 
to  look  for  a generally  valid  answer.  The  problem  wi 
Investigated  for  the  simple  case  of  only  three  simultaneous  reflections 
and  a limited  number  of  g.b  combinations.  In  addition,  one  reflection 
Ewald  sphere  (s  > 0).  Since  a strong 
ier  reflections  can  in  practice  only  be 
of  hexagonal  symmetry  (rectangular 
im  cases),  two  220  reflections  of  the  111 

> deviation  p. 


shown  in  Figures  3d 


realised  in  diffraction  patterr 

diffraction  plane  were  chosen  a 
forming  an  equiaxed  triangle  together 

Calculations  could  be  carried  out  up  to 

Runge-Kucta  integration.  Some  of  the  results  a: 
through  36  which  represent  the  conditions  g'b  - 1,1,  g*b  ■ 1,-1  and 
g-b  - 1,2,  respectively  (the  case  of  g.S  - 1,1  will,  for  three  strong 
beams,  in  practice  lead  to  a six-beam  orientation  due  to  Che  symmetry 
of  the  111  diffraction  pattern).  In  all  three  figures,  bright  field 
(a)  and  dark  field  (b)  profiles  are  calculated  for  six  different  values 

1,  -1  and  2,  respectively,  in  the  three  figures)  ranging  from  w-  ■ 0 

for  the  Runge-Kutca  routine  to  complete  the  computation  of  one  profile 


profile  points  indicated  no  major  changes  w: 


Three-beam  intensity  profiles  across  a screw  disloca- 
tion for  a foil  orientation  indicated  by  the  diffrac- 
tion spots,  g ■ 022  (g-S  - 1)  always  at  the  reflect- 
ing position,  while  w is  varied  for  g “ 202  (g.'i>  - 1), 
Bright  field  in  (a),  dark  field  in  (b).  t » 3.0  £ , 


(a) 


As  may  have  been  expected,  the  series  of  profiles 
34,  35  and  36  look  quite  complicated,  particularly  when  be 
beams  are  excited  to  a comparable  degree  (profiles  in  (a) 
toward  the  upper  left).  A preference  for  three-beam  doub! 
a particular  combination  of  g.S  products  (see  Howie  and  Wh 
cannot  be  concluded  from  this  limited  evidence. 


i diffracted 


quite  fascinating 


figures  at  larger  w for  the 
in  a different  light  when  tt 


■n  profiles  in  all  three 
with  a general  decline  of  the 

e last  profiles  in  Figures  34  through  36  are 
brighr-fleld  calculations,  for  example, 


show  that  the  third-beam  influence,  even  at  w„  ■ 4.0,  is  still  strong 
enough  to  produce  a profile  which  essentially  represents  the  opposite  ol 
the  two-beam  Image.  It  Is  then  quite  surprising  that  the  Intensity 
i,  being  subject  to  such  a profound  third-beam  effect. 


Without  further  study,  the  only  definite  conclusion  which  c. 
drawn  from  the  limited  number  of  calculations  presented  above  reduci 
the  fact  that  a possible  third-beam  influence  cannot  be  Ignored  as  : 
as  w3  is  smaller  than  or  equal  to  5.0.  In  some  cases  (smaller  atom: 
spadngs  in  heavier  metals,  or  lower-order  reflections) , a value  of 


le  reciprocal  let  rice 
■ remembered  that  all 

le  diffraction  spot 


v - 5.0  may  represent  a very  large  deviation  of  t 
point  from  the  Ewald  sphere.  It  must  therefore  l 
cany -be- am  profiles  apply  only  to  the  particular  c 
reflections  on  which  the  calculation  has  been  bas 
sentative  parameter  may  be  the  distance  between  t 
and  the  Kikuchi  line  for  the  third  beam  which  in  the  case  of  w > ! 
the  profiles  of  Figures  34-36  would  be  equal  to  about  one  fourth  < 
220  reciprocal  lattice  vector  (this  may  be  compared  to  a distance 
one  reciprocal  lattice  vector  for  220  in  the  systematic  case  with 

increasing  foil  thickness,  as  discussed  in  Section  4.4.1.  should  l 


Although  dynamical  contrast  calculations  often  show  the  danger 
involved  in  projecting  results  obtained  for  one  set  of  conditions  to 
hold  under  a different  sec,  it  may  be  speculated  for  a moment  chat  three- 
beam  profiles  would  remain  Independent  of  the  g -b  product  down  to  the 
region  of  systematic  beam  deviations  (w  19.1  for  the  220  row  in  aluminum). 
In  this  case,  many-beam  reflections  might  become  independent  of  the 
direction  of  diffracting  vectors,  and  many-beam  corrections  to  two-beam 
images  may  only  be  a function  of  the  number  of  beams  and  possibly  the 
order  of  reflections.  Also,  such  corrections  could  be  expected  to  hold 
for  quite  arbitrary  dislocation  types.  Though  nothing  but  a speculation 

stimulate  further  study  of  three-beam  contrast. 

Besides  the  purely  theoretical  calculations,  matching  was  tried 
for  profiles  generated  by  three  or  more  strong  beams.  Figure  37b 


on  Intensity  profile  ai 


luo-beam  profile  shown  for  comparison, 
intensity  scale  0.1  per  division. 


Calculated 


113 


represents  the  Intensity  profile  across  the  screw  dislocation  shown  in 
Figure  37a,  produced  by  essentially  three  strong  beams.  The  match 
with  the  calculated  profile  in  Figure  37c  is  quite  good,  and  the  two- 
beam  profile  added  for  comparison  in  Figure  37c  indicates  the  strong 
effect  of  the  weaker  beam  in  this  case.  In  other  situations,  however, 
agreement  between  computed  and  experimental  profiles  was  very  poor, 
despite  a careful  evaluation  of  deviation  parameters  using  the  orien- 
tation program  developed  in  Appendix  A2.  Although  many-beam  matching 
is  not  believed  to  pose  a serious  difficulty  in  principle,  the  greater 

values  seers  to  require  a much  larger  amount  of  calculation  than  is 
necessary  for  a two-beam  case. 


The  image  scanning  system  built  into  the  electron  microscope 
°ffers  the  advantage  of  sampling  the  electron  intensity  directly, 
without  relying  on  a secondary  process  of  intensity  conversion  which 
the  microdensitometer  requires  in  form  of  the  photographic  place.  In 
Figure  38,  intensity  profiles  across  a stacking-fault  dark  field  image 
in  Cu-Al  are  shown,  produced  with  the  microdensitometer  (Figure  38a)  as 
well  as  by  direct  intensity  scanning  (Figure  38b).  The  two  micro- 
densitometer traces  in  Figure  38a  represent  two  different  photographic 
plates  which  were  exposed  for  30  sec  (solid  line)  and  8 sec  (dotted 
line)  in  the  electron  microscope.  Similar  differences  in  traces  of  an 
identical  image  were  obtained  with  constant  exposure  but  different 
developing  times.  Each  trace  would  lead  to  a somewhat  different  set  of 


responsible  ft 


values  for  relative  peak  height  and  half-peak  width,  so  that  an  accurate 
calibration  of  the  photographic  process  will  be  necessary  when  using 
the  microdensitometer  for  image  analysis.  Although  the  Philips  EM  200 

to  determine  the  correct  exposure  time,  practical  difficulties  may 
arise  when  most  of  the  electron  Intensity  is  concentrated  in  an  area 
which  is  small  in  comparison  to  the  screen  site.  Such  a case  is  often 
ark  field  images  at  lower  magnifications. 

‘ther  hand,  direct  intensity  scanning  suffers  from  the 

the  large  noise  in  the  signal  current  (see  Figure  38b). 
The  noise  level  is  strictly  a function  of  total  electron  intensity  in 
the  image  which  at  present  requires  the  scanning  system  to  operate  at 
either  maximum  beam  divergence,  relatively  low  magnification,  or  in 
rather  thin  specimens.  The  present  system  is  therefore  clearly 
inferior  to  the  microdensitometer  in  high-resolution  work. 

The  problem  of  increasing  the  total  electron  intensity  for  a 
given  image  point  during  scanning  has  been  attacked  in  two  ways.  First, 
the  image  movement  across  the  scanning  aperture  was  made  discontinuous, 
and  Che  number  of  electrons  which  passed  through  the  aperture  in  a given 
time  (of  the  order  of  some  seconds)  at  each  calibrated  increment  was 
recorded  by  a digital  counter.  Figure  39  shows  an  example  of  this 
approach,  where  counts  were  taken  across  a dissociated  dislocation  in 

image  peaks  (the  large  fluctuations  of  the  continuous  scan  are  typical 
for  a low-intensity  image  in  a relatively  thick  foil).  In  general. 


Intensity  profiles  of  an  extended  dislocation 
in  a Cu-13.3  at  1 A1  alloy.  Solid  line  represents 
direct  scan  across  the  low-intensity  image.  Data 
along  broken  curve  originate  from  discontinuous 

the  image.  Accumulated  counts  between  1000  and 
5000  per  image  point. 


however,  this  method  proved  to  be  impractical  due  Co  the  image  instability 
during  the  necessarily  long  period  of  time  to  complete  the  profile. 

The  second  approach  utilized  the  variable  frequency  generator 
of  a signal  averager  whose  saw-tooth  profile  caused  the  image  to 
oscillate  back  and  forth  across  the  scanning  aperture.  During  each 
sweep  period,  the  number  of  signal  electrons  per  time  interval  (250 

counted  and  stored,  and  the  accumulated  counts,  as  a function  of  distance, 
at  any  time  displayed  on  a cathode  ray  tube.  By  this  method  of  signal 
averaging,  image  peaks  were  expected  to  sharpen  with  time,  due  to  the 
relative  reduction  in  random  noise  after  repeated  sweeps.  Sweep 
periods  could  be  chosen  over  a range  from  one  second  to  a few  micro- 

A possible  reason  for  this  failure  may  have  been  a memory 
effect  (hysteresis)  in  the  scanning  colls.  Such  an  effect  was  found 
when  an  ordinary  profile  was  repeated  in  the  backward  direction  although, 
within  the  accuracy  of  the  intensity  measurement,  image  peaks  would 
appear  at  identical  positions  when  scans  were  repeated  in  the  same 

registered  while  the  image  was  deflected  in  one  direction,  due  to  the 


(30  sec)  to  complete 
sampled.  From  this  p 


ordinary  scan. 


accuracy  to  analyze  the  Images  of  steeply  Inclined  dislocations  dis- 


cussed in  Section  4.2.  The  chsnneltron  scanning  system,  however, 
offers  the  advantage  that  a magnified  image  of  the  defect  (about  8X 
with  respect  to  the  photographic  plate)  can  be  observed  as  it  crosses 


the  aperture,  whereas  the  microdcnsitometer  requires  alignment  of  the 
photographic  plate  at  a rather  large  distance  from  the  light  slit.  The 
present  limitation  of  the  channeltron  detection  device  is  solely  due 
to  a large  hole  in  the  viewing  screen  above  the  scanning  aperture  (see 
Figure  5),  yet  the  correlation  between  scanning  position  and  defect 
image  was  still  more  accurate  than  with  the  microdensitometer. 

A final  criterion  of  practical  importance,  which  favors  the 
direct  intensity  scanning  method,  concerns  the  time  factor  involved  in 
electron  microscopy  image  interpretation.  The  microdensitometer  can 
only  furnish  results  after  the  lengthy  procedure  of  photographic 
processing,  at  which  time  the  original  image  in  many  cases  is  no  longer 
available,  although  verification  of  a particular  feature  or  additional 
analysis  may  be  desirable.  Clearly,  the  Instant  date  supply  from  the 
direct  scanning  system  provides  such  an  opportunity. 

Thus,  it  may  be  concluded  from  the  experimental  evidence 

represents  the  better  choice  for  generating  image  profiles  if  the 
analysis  con  be  carried  out  at  relatively  low  magnification  (<30000  X 
on  the  photographic  place)  and  only  moderate  resolution  requirements. 
Without  considerable  improvement  of  the  signal-to-nolse  ratio  in  the 

both  x-scale  end  intensity  scale,  should  at  present  employ  the  micro- 


CHAPTER 


CONCLUSIONS 


A first  step  has  been  made  toward  a more  quantitative  evalua- 

experimental  dislocation  profiles  have  been  analyzed  to  determine  the 
sign  of  the  Burgers  vector,  and  to  distinguish  between  the  effect  of  a 
traction-free  or  a rigid  foil  surface  on  the  displacement  field. 

Such  an  approach  illustrates  the  method  which  is  required  when 
an  unknown  defect  must  be  identified  by  selecting  the  appropriate  form 
of  its  displacement  field.  The  calculations  in  Section  4.2  have  shown 
that,  in  general,  one  will  have  to  consider  a rather  large  number  of 

of  contrast  parameters.  It  is  therefore  desirable,  for  matching  pur- 
poses, to  rely  on  as  much  experimental  information  from  Che  image  as 
possible  because,  in  profile  calculations,  any  one  feature  detsll  may 
possibly  be  generated  from  different  displacement  fields  at  different 

1 profile  scanning.  In  the  case  of  dislocations,  it  eliminates 
accurate  measurement  of  the  depth  position  in  the  foil. 


Secondly,  it  combines  the  total  image  information  in  the  form  of  a 
single  picture.  In  principle,  che  total  information  can  also  be  produced 
by  a number  of  profile  scans  taken  at  various  positions  along  the  defect. 
However,  beside  the  time  limit  imposed  on  scanning  by  image  instabilities, 
the  practical  difficulties  should  be  emphasized  which  Che  human  mind 


when  constructing  Che  total 
individual  intensity  profiles. 


reasonably  well  by  the  two-dimensional  image  simulation,  some  doubts  may 
be  left  about  the  correct  form  of  the  displacement  field,  which  can 
only  be  removed  from  a knowledge  of  the  quantitative  intensity  distri- 
bution at  some  characteristic  positions  along  the  image  (e.g.,  consider 
the  oscillating  dislocation  concrast  in  Section  4.2).  It  is  therefore 
suggested  that  the  interpretation  of  unknown  defects  from  strain  contrast 
should  initially  proceed  by  using  the  image  simulation  method,  while 
quantitative  profile  scanning  would  serve  as  refinement  to  this  technique. 
The  discussion  in  Section  4.5  should  be  consulted  in  order  to  choose  the 
optimal  procedure  for  generating  intensity  profiles. 

Quantitative  contrast  evaluation  as  described  in  Section  4.3 
has  revealed  a tendency  for  image  broadening  during  standard  operation 
of  the  electron  microscope.  Although  the  interpretation  of  unknown 
defects  could  conceivably  ignore  absolute  profile  widths,  image  broaden- 
ing would  play  an  important  role  in  quantitative  measurements  of  coherency 
strains.  In  the  analysis  of  Ashby  and  Brown/®  these  strains  are 
directly  proportional  to  the  image  width  of  the  defect.  The  study  of 
dislocation  width  should  therefore  be  continued  in  order  to  obtain  more 
accurate  estimates  of  the  extent  of  image  broadening  as  a function  of  the 
operating  reflection.  This  should  include  a more  extensive  matching 

if  beam  divergence  and  Inelastic  scattering  on 
sis  to  provide  an  interesting  study  for  further 


image  broadening  s. 


analysis  of  the  mechanism  of  strain  contrast.  Clearly,  electron  energy 
analysis  promises  considerable  advantages  with  respect  to  the  inelastic 
scattering  effects  in  such  a study,  and  has  in  fact  received  widespread 


The  accuracy  of  any  evaluation  of  strain  contrast  will  depend 
on  effective  values  of  the  anomalous  absorption  parameter.  Zt  is 
therefore  suggested  that  the  study  of  absorption  from  oscillating 

extended  to  a more  quantitative  investigation  by  generating  intensity 
profiles  across  the  black-white  dislocation  images.  For  this  purpose, 
the  viewing  screen  above  the  electron  detection  system  in  the  micro- 
scope should  be  redesigned  in  order  to  allow  better  monitoring  of  the 
precise  scanning  position.  If  the  results  suggested  in  Figure  18  can 
be  confirmed  quantitatively,  accurate  absorption  parameters  may  become 
available  for  strain  contrast  experiments  by  simply  recording  intensity 
profiles  across  thickness  extinction  fringes  which  can  be  found  in  meet 
foils  produced  by  electrolytic  chinning  methods. 

The  effect  of  foil  thickness  on  the  magnitude  of  the  effective 
anomalous  absorption  parameter  may  prove  useful  for  a better  understanding 
of  Che  contribution  of  inelastic  scattering  to  strain  contrast.  Such 
a study  would  therefore,  again,  benefit  from  the  presence  of  an  energy- 
analysing  section  in  the  microscope.  Different  absorption  parameters 
have  in  this  way  been  measured  by  Takagl  and  Ishida^^  from  "no-loss" 
electron  images  of  thickness  fringes. 

Mony-beam  effects  reported  in  Section  4.4  have  been  surprisingly 
strong,  both  with  respect  to  systematic  and  nonsyscematic  reflections. 


As  mentioned  previously,  many-beam  calculations  have  become  a popular 
research  activity  due  to  the  technological  advances  of  high-voltage 
electron  microscopy.  But  even  if  the  two-beam  orientation  is  retained 
as  basis  for  experimental  strain  contrast  analysis  at  100  kV,  these 
first  results  suggest  chat  the  theoretical  investigation  of  weak-beam 
effects  be  continued,  along  the  lines  discussed  in  Section  4.4,  in 
order  to  determine  the  reliability  of  two-beam  calculations  in  particular 
contrast  situations.  For  this  purpose,  the  numerical  technique  of 
integrating  the  intensity  equations  (2.30)  in  Section  2.2  should  be 
reexamined . 


major  roads  on  which  electron  microscopy  should  proceed  in  the  future 
in  order  to  extend  its  usefulness  for  applied  science  and  engineering. 
The  emphasis  on  quantitative  results  in  this  thesis  represents  a new 
development  in  the  evaluation  of  strain  contrast.  Beside  proving  the 
feasibility  of  such  an  approach,  its  promise  of  supplying  useful 
information  has  been  demonstrated  in  some  particular  contrast  applice- 


ic  can  be  developed  as  a new  tool  for  electron  microscopy. 

The  results  of  this  thesis  may  be  summarised  as  follows: 

(1)  The  experimental  accuracy  of  measured  contrast  parameters 
in  general  will  not  be  sufficient  to  determine  uniquely 
the  shape  of  a dislocation  intensity  profile.  Extensive 
contrast  calculations  therefore  constitute  a prerequisite 
for  the  interpretation  of  unknown  displacement  fields. 


(2)  Quantitative  profile  analysis  c. 


of  image  broadening  seems  to  be  a function  ol 
reflection. 

(4)  the  relative  position  of  a dislocation  image 
to  the  true  dislocation  line  position  agrees 

(5)  The  effective  anomalous  absorption  parameter 
of  foil  thickness,  and  seems  to  agree 
from  the  intensity  distribution  across  thickness  extinction 

(6)  Strong  contrast  effects  from  weak  systematic  and  non- 
systemaeie  reflections  must  be  expected,  at  least  in  thin 

locations  in  HCP  as  well,  and  allows  in  situ  analysis 
of  the  dislocation  character  on  the  viewing  screen  when 
slip  plane  and  foil  plane  coincide. 
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APPENDIX  2 

CALCULATION  OF  DEVIATION  PARAMETERS  FOR  NON- 
SVSTEMATIC  REFLECTIONS 


The  basic  idea  for  Che  following  derivation  is  to  describe  the 
oriencation  of  the  incident  beam  with  respect  to  the  foil  by  the  posi- 

the  reciprocal  lattice  plane  with  the  Ewald  sphere  (Figure  40).  It 

circle,  does  not  lie  in  the  plane  of  the  diffraction  pattern  except 
for  the  case  of  the  so-called  symmetrical  orientation  (beam  direction 
parallel  to  Che  zone  axis),  when  it  degenerates  to  a point. 

For  the  simple  case  of  systematic  reflections  only,  the  problem 
con  be  solved  in  two  dimensions  (Figure  41) , where  the  Ewald  circle 


P„.  Bragg's  law  is  satisfied: 


and  a distance,  defined  as  that  by  which  such  a point  is  separated  from 
the  plane  C*  parallel  to  the  reciprocal  lattice  plane,  can  be  expressed 


Bragg  angles 


lattice  plane,  uw,  with  the  Ewald  sphere. 


110 


diffraction  pattern.  All  reflections 
from  C have  identical  deviation  parameters. 


o points  which  mi 


Inserting  (A2.8) 


expressions 
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s plane  of  diffraction,  the 


e applied  to  determine  the  ci 


now  be  obtained  by  noting 


reciprocal  lattice 


Thus,  the  point  C is  finally  obtained  with  Indices 
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101  FORMAT ( 1 2 > * 

111  FORMAT t F10.A) 

201  FORMAT [‘OCENTER  0 
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27  CONTINUE 

NPLOT =BPLUT ( I ) 

SHIN=0.0 

SMAX=0.5 

CALL  PLOT2  ( GRAPH, XMAX, XM IN, SMAX, SM IN ) 

2016  FORMAT  1 1111 ) 

CALL  FPLUTA  (1,1H$) 

31  CONTINUE 
7 CONTINUE 

501  WRITE  (6,2511 

502  STOP 


NUMERICAL  INTEGRATION  OP  DIFFERENTIAL  EQUATIQNS- 
FOURTH  OROER  RUNGE-KUTTA  METHOD 
SUBROUTINE  RUKU(H,M,Y,0,E,N,J) 

DIMENSION  YI300), 0(501 
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